Cluster states of Fermions in the single /-shell model 



in 
o 
o 

> 
O 

o 



Andras Csordafl 

Research Group for Statistical Physics of the Hungarian Academy of Sciences, 
Pdzmdny P. Setdny 1/A, H-1117 Budapest, Hungary 

Eva Sz6k(li| 

Department of Physics of Complex Systems, Eotvos University, 
Pdzmdny P. Setdny 1/A, H-1117 Budapest, Hungar'^ 

Peter SzepfalusjH 

Department of Physics of Complex Systems, Eotvos University, 
Pdzmdny P. Setdny 1/A, H-1117 Budapest, Hungar-^ and 
Research Institute for Solid State Physics and Optics, P. O. Box 49, H-1525 Budapest, Hungary 

(Dated: February 2, 2008) 

The paper concerns the ground state structure of the partly filled i-shell of a fermionic gas of 
atoms of spin s in a spherically symmetric spin independent trap potential. At particle numbers 

= n(2s + 1), n = 1, 2, . . . , 2/ + 1 the basic building blocks are clusters consisting of (2s + 1) atoms, 
whose wave functions are completely symmetric and antisymmetric in space and spin variables, 
respectively. The creation operator of a cluster is constructed whose repeated application to the 
vacuum leads to the multi-cluster state. Ground state energy expressions are derived for the n- 
cluster states at different I, s values and interpreted in simple terms. 
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I. INTRODUCTION 

The many body problem as applied to finite systems 
has a long history in atomic and nuclear physics [ll Q, H, 
Q. One of the central problems has been the nature of 
the ground state in case of a partially filled shell. 

When we consider atoms (a gas at zero temperature) 
in an external potential the possible behaviors are quite 
rich. One can assume that the collision between atoms 
does not excite internal degrees of freedom and the atoms 
can be regarded as structureless objects whose spins are 
fixed being in a definite hyperfine state. The external 
potential can be supplied by a magnetic or optical trap 
[5|- In the past few years a very intensive research has 
been continued in case of Bose-particles both theoreti- 
cally and experimentally in a variety of such systems. 
Fermion systems along these lines have been less stud- 
ied until recently, but important achievements have been 
already available and one can be sure that rapid devel- 
opment will continue in the future. In particular the 
achievement and study of the superfluid state has become 
one of the frontiers in physics (see for Reviews |& iTj, |8|, i^] , 
which contain references to earlier works). 

In this paper we treat fermionic systems. The exter- 
nal potential will be assumed to be spherically symmet- 
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ric and the interaction between atoms will be described 
by a spin independent 5 function like attractive pseu- 
dopotential. Within the single ^-shell model of fermionic 
atoms with arbitrary large spins we study the ground 
state properties of the gas. 

We concentrate on cluster states containing particles 
of numbers Nn = n{2s -\- 1) where s is the spin of the 
particles and n = l,2,...,2Z-|-l. The existence of such 
clusters have been pointed out by us previously [Tol |. In 
the present paper we extend the investigations in several 
directions. It is advantageous to use the second quan- 
tized representation, which makes possible, among other 
things, to obtain the multi-cluster states by repeated ap- 
plications of of the cluster creation operator to the vac- 
uum. Besides numerical calculations extensive analytical 
studies are carried out. It is shown that to a good ap- 
proximation the hamilton operator can be replaced by a 
model one built of the the operator of the particle num- 
ber, of the quadratic Casimir operators of the groups 
SU{21 + 1) and SO{2l + 1). It is shown that within 
the model the clusters can be conceived as interacting 
(Cooper) pairs if s > 1/2, the interaction being of statis- 
tical origin. The model Hamiltonian coincides with the 
true one for / — 1,2 (The special cases mainly investi- 
gated in [loj). 

For the model hamiltonian the cluster states are ex- 
plicitly given and it is pointed out that the energy of the 
ri-cluster state can be written as the sum of energies of 
the clusters proportional to n and an "interaction term" 
term between the clusters proportional to n{n — l)/2. It 
is found that the second term disappears when s — 1/2. 
For s = 1/2 it was shown by Racah in 1952 [ij that the 
ground state consists of independent pairs. 
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Investigations of particles in an open shell has been an 
important area in atom and nuclear physics. As discussed 
above in the trapped gas of Fermi atoms new features ap- 
pear due to the fact that the spin of the particles can be 
higher than 1/2. Comparing with the situation in the 
atom shell a further important difference is that instead 
of the long ranged Coulomb force between the electrons 
the atom-atom interaction which has to be considered is 
short ranged, while comparing the situation with that of 
an open neutron shell in a nucleus an important differ- 
ence is that in our case the open shell is an l-one (since 
no-spin orbit interaction is present) as contrasted to the 
j-one in the nucleus. 

Though our main purpose in this paper is to enlarge 
the picture we have about the dynamics of fermionic par- 
ticles in a partly filled shell, a few words about the rel- 
evance of the model for physically realizable situations 
in case of trapped gases are in order. Obviously optical 
traps are the suitable ones which allow the free rotation of 
the spins. It is assumed that the external potential con- 
tains besides the confining potential (which is typically 
of a harmonic oscillator type) the mean field of the atoms 
building the closed shells. There are two conditions then 
to be fulfilled in order that the single-l shell model apply. 
Firstly, the mean field due to the closed shells should be 
strong enough that the possible degeneracy of the lev- 
els of different I values be lifted considerable. Secondly, 
the characteristic interaction energy of the atoms in the 
partly occupied shell should be smaller than the relevant 
level distance in the external potential. Both require- 
ments can be satisfied simultaneously if the strength of 
the interaction between the atoms is weak and the num- 
ber of the atoms in the trapped gas is sufficiently large 

Elll. 

The paper is organized as follows. In Section |n] 
the Hamiltonian is written in terms of irreducible ten- 
sor operators and a model Hamiltonian is introduced 
built from the operator of particle number and from the 
quadratic Casimir operators of the groups SU{21 + 1) and 
SO {21 -1-1). The Section Hill is devoted to the extensive 
investigation of the one-cluster state. The symmetrizer 
playing an important role in creating the cluster states 
is analyzed in detail in Section IIVI Energy eigenvalues 
of the model Hamiltonian for cluster states generated by 
the cluster creation operator are calculated in SectionM 
It is shown in Section IVII that the expectation value of 
the Hamiltonian agrees with that of the model one for 
cluster states. Section IVIII contains the summary and a 
discussion of the results. Appendices ^ and ^ present 
certain steps of proofs outlined in the main text. Ap- 
pendix El discusses the special case when the spin of the 
particles is 1/2. 



II. FORMULATION 

Let us assume that the trap potential is spherically 
symmetric. At zero temperature the open shell consists 



of fermions of spin s which partially fill the (n, I) shell. 
The one particle normalized wave functions are given by 

Here is the discrete spin variable and xti"^) the nor- 
malized spin eigenfunction. In the model the quantum 
numbers (n, /, s) are fixed, m and v can take the values 
m = —I, — I + 1, ... J and v = — s, —s + 1, . . . ,s and the 
particle number N can vary between and (2/-l-l)(2s+l). 
The functions £3) are eigenfunctions of the one-particle 
Hamiltonian, which contains the kinetic energy term, the 
trap potential and the averaged field potential of the 
closed shells. The interaction between the particles in 
the partially filed shell is given by 

A ^ 

H,^t = --Y, Sir,-v,), A>0 (2) 

i,j=i 

corresponding to a spin independent s-wave scattering 
with negative scattering length. Our aim is to diagonalize 
(0) on the fixed basis ((TJ. 

Let us denote the operator which annihilates a parti- 
cle with quantum numbers (n, /,m, s,i^) by a^.u- The 
Hamiltonian in second quantization reads as 

Hint = hEo/TT — — ^ ^ ^ ^ /mi,m2;m3,m4 

^li^2 1/1,1/2 

where Eq is the characteristic energy 

Eo = \ \RnAr)\^r^ dr, (4) 
Jq 

and / can be expressed in terms of the Wigner-Sj symbols 
|14| in two equivalent forms 

fm,,n.r,^,,m, ^ J df7 < (f7)i;[:^ (f^)^^ (f^X (f7) 



X 




For notational simplicity we introduced the symbol [. . .] 
defined by 

H = (2p+1). 
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The dimensionless Hamiltonian h (see Eq. Q) can be 
written as 



8 ^ \^ / 



where 



Bl= i--^)'-"B,,MBL, 



-M 



(8) 



M=~L 



(in eq. ([TJ only the terms L even remain). The operators 
13 L M defined as 



I s 
m= — l iy= — s 

I I L 

m M — m —M 



X 



^ a+^am-Af,y (9) 



are spin scalars and irreducible tensoroperators with re- 
spect to angular momentum IJ. Special cases are L = 



N 



0,0 



(10) 



and L = 1: 



V3 



Si,±i = T 



vw+m 

V3 



(11) 



The operator _Bo,o commutes with all the others and the 
operators B^^m for L > 1 form a Lie-group, which is 
isomorph to SU{21 + 1) with the commutators 



E VMlW] 



L".M" 



1 - (-1)^+^'+^" (-1)^^" 



L L' L" 
M M' -M" 



L L' L" 
I I I 



Bl",m" 



(12) 



({. . .} denotes the Wigner-6j symbol). Due to the special 
form of the structure coefficients the operators -Bl.m for 
odd L form a subgroup, which is isomorph to 50(2/ -I- 1). 
This latter also has a subgroup 50(3) 1^ spanned by 
Bum, M = 0,±1. The Casimir-operator of SU{21 + 1) 
is 



21 



(13) 



and that of 50(2/ -I- 1) 



21-1 



a = E ^ 



(14) 



L = l 
L:odd 



The hamiltonian h defined in Eq. Q commutes with , 
L^, 5^, Sz and with the Casimir operator Cu- In fact, 
the Wigner-3j symbol vanishes for L odd, thus in Eq. {Tj) 
the sum over L runs over even values of L. In the special 
case / = 1, because the operator — Cu — Co, ^"^d 
case / = 2, because of the accidental coincidence 
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h can be expressed entirely in terms of Cu and Co- 
For / > 2 this is not true anymore, and furthermore Co 
does not commute with h. 

However, let us consider the model-Hamiltonian 
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hrr,. = ^-^N - — 



[I] 



8 4(2/ + 3) 



(Cu - Co) , (15) 



This model hamiltonian has the following important 
properties: In the special cases 1 — 1,2 the two Hamil- 
tonians h and hm agree. For / > 2 the two hamiltonians 
are different and we write for the difference 



R = h — hr, 



E 



/ / L 




(2/ + 3)[/] 



Bi 



L=l 



(16) 

In the next sections we show that the averaged value of 
R with the ground states of hm is zero if the particle 
number equal to 

Nn = n{2s + 1), n = 0,l,...,(2/-f 1). (17) 

i.e., N is a multiple of (2s -1-1). In ^3 have reported 
that a clusterization in the wave function occur at par- 
ticle numbers Nn and that the ground state energy per 
particle has local minima at the n-cluster particle num- 
ber Nn- We also show that h and hm share the same 
ground states for even particle numbers if s — 1/2. 



III. CLUSTER STATES 

Making a full numerical diagonalization of the hamil- 
tonian (jZj) is not easy. Fock- vectors have (2/ -I- l)(2s + 1) 
slots, and in each slot there is a zero or 1 due to the 
fermionic character of the problem. The dimension of 
the full Hilbert-space is 2(2'+i)(2'*+i). The spectra do not 
depend on Lz and Sz, thus we can restrict ourshelves 
to the fermionic sectors Lz — 0, Sz = (even particle 
numbers) or Lz = 0, Sz = 1/2 (odd particle numbers). 
Conserved operators such as 5^, L^, N and Cu makes 
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the numerical problem block-diagonal, but still the com- 
puter time and storage required grows exponentially fast 
as soon as we increase the open shell quantum numbers s 
or I. This motivates our analytical approach besides the 
numerical efforts. 

In we have investigated the ground state wave 
function at the particle number A'^i in first quantization. 
In second quantization the corresponding wave function 
reads as 



I Id) 



k 



CkSn 




|0), 



k = (mi, . . . ,m2s+i). 



(18) 



with some coefficients Cfc. Here the symbol Sm is an 
operator which symmetrizes its argument with respect 
to the indices mi, ... , m2s+i (see Sec. II V|) . 

For notational simplicity let us introduce the integer a 
by 



{2s + 1) 

The ground state of hm, explicitly given by 



(19) 



(20) 



(see Sec. 0. The operator Qqq in (pn|l shall play a 
central role in the following analyzis. It creates a pair 
state from the vacuum with quantum numbers L = S = 
with fixed open shell quantum numbers {I, s): 



1 



io.o — 



Special cases of (|TH|l are 

\ld) = \ld)o, I ' " 



-u+t-m + + 

m.i' —m. — jy 



1,2, 
1/2, 



s: arbitrary 
^: arbitrary 



(21) 



(22) 



since h and h„i coincides for / — 1,2. The equality for 
s = 1/2 will be considered in Appendix lO 

The first example when (|22(l is not valid occur for I = 
3, s = 3/2. = A^i = 4 is the first-cluster particle 
number. The ground state is in the L = Lz = S — Sz— 
fermionic sector. In this sector there are five orthonormal 
basis vectors. The matrix elements of h on this basis are 
the following 



/O 














_ 14 

11 







4886 
' 2145 



35 
6 



The true ground state energy for = 4 comes from the 
lowest 2x2 block-diagonal and has the value 



E 



7(1657 -hV537729) 
2860 



-5.85038, 





s=l/2 


s=3/2 


s=5/2 


s=7/2 


1=1 


-0.75 


-3.30 


-7.65 


-13.80 


1=2 


-1.25 


-4.64 


-10.18 


-17.86 


1=3 


-1.75 


-5.85 


-12.32 


-21.17 


1=4 


-2.25 


-6.99 


-14.26 


-24.12 


1=5 


-2.75 


-8.09 


-16.09 


— 


1=6 


-3.25 


-9.17 


-17.83 


— 


1=1 


-0.75 


-3.30 


-7.65 


-13.80 


1=2 


-1.25 


-4.64 


-10.18 


-17.86 


1=3 


-1.75 


-5.83 


-12.25 


-21.00 


1=4 


-2.25 


-6.95 


-14.11 


-23.72 


1=5 


-2.75 


-8.04 


-15.86 




1=6 


-3.25 


-9.10 


-17.55 





TABLE I: Ground state energies of h (upper part) and hm 
(lower part) up to two decimal digits precision at the one- 
cluster particle numbers Ai = (2s + 1). 



which is quite close to ^,5,5 — —35/6 « —5.83333. Nu- 
merically hi,j in the same 2x2 block is 



-2.27786 -0.246813 
-0.246813 -5.83333 



and the lowest energy state belongs to the eigenvector 

0.0690865 ^ 1-23) 

It is interesting to present the matrix elements of C„ 

264 



C{SUi7)),., 



and that of Co 



C(50(7))„- = 





264 
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There is a small admixture of two vectors belonging to 
different eigenvalues of C'o and the dominant vector-part 
of the eigenvector belong to such a vector, which is an 
eigenvector of Co with eigenvalue zero. Analyzing further 
this fifth vector it turns out that it is still given by H20|) 
with ^ = 3, s = 3/2. The average value of h by this vector 
is equal to that of hm- Similar property will be proven in 
Sec. IVII also for multi-cluster states. It is important to 
stress that the fourth vector and the true ground state 
Vg are also a symmetrized states with S = 0. 

In Tabled we show the numerically calculated ground 
state energies oi h ^ and h„i H15|l . From the data it 
is clearly seen that for Z = 1,2 or for 5 = 1/2 Eq. H22|) 
exact and for all the other cases \ lcl) |lcZ)o is a rather 
good approximation. This shows the importance of the 
model hamiltonian hm- The dominating contribution to 
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the energy (lower part of Table HJ can be written in the 
form 



(2) (^i<^ - 1) r 

ei = -ae^''' — -S, 



where e'^' and S are independent of the spin. Their ex- 
plicit expressions will be derived in Sec. The first 
term gives the energy of a independent pairs while the 
second term lowers this energy, indicating clearly that the 
cluster wave function gives lower energy than the wave 
function of independent pairs. Note that an eigenfunc- 
tion of hm exists leading to the eigenvalue — CTe*^^-* if the 
inequality a < 21 + 1 is fulfilled. For spin one half parti- 
cles one has only of course, the first term, the energy of 
a single pair. 

Both \lcl) and |lcZ)o involve the operator Sm, there- 
fore, let us study the properties of the symmetrizer Sm- 



IV. THE PROPERTIES OF THE 
SYMMETRIZER Sm 

Let us introduce the symmetrizer symbol Sm by 
the following properties: i) Sm is linear for its argu- 
ment, ii) Sm symmetrizes any pure operator-product 
^mui^i ■ of creation operators i = l,...,p 

with respect to all without changing the order of spin 
projections (z^i, . . . , Vp) including the combinatorical nor- 
malization. As an example: 

'S'm(a^3y2'^l^l/2"o^-l/2) ~ 3 ('^^3/2'^^ 1/2^^(1 -1/2 
+ '^l,3/2'^0,l/2'^0,-l/2 + '^0,3/2'^0,l/2'^l,-l/2) ■ 

{After symmetrization, of course, one can change the 
order of creation operators as we wish using the anti- 
commutativity of creation operators) . 

Very important special case is p — 2s 4- 1. Writ- 
ing down all the terms in Sm{a+^^ ,,^ ■ ■ ■ a^a.+i^i^s^+i) 
turns out that it must be antisymmetric in all spin index 
(j/i, . . . . i^2s+i)- Therefore, by inspection we have 



Smio-mi,ui ■ ■ ■ '^m23 + i,i^2a + i) ^ ^^i^i ,. . .,!^2s+l 
xSm{0'mi,s ■ ■ ■ ^m2s + i,-s)' 



(24) 



where ei/-^^...^u23-i-i is the antisymmetric tensor with the 
convention es,s-i,...,-s — 1- Further expansion is possible 
for the combination 



Smio^mi, 



. a. 



+ 1)! ^ 



m2.+ i,-s/ (2S+1)! ^ 

Vl 1^23 + 1 



(25) 



■^i/i ,. . . ,i.'2s-fi mi ,1/1 ■ • ■ '^m2s4.i,f2s + i 

Next, we enumerate some properties of useful in the 



following. If 

■^a ^ ^ 9rti!n^m,u^n,~i'^ Ck. = 1, . . . 

(26) 



(a) = (a) 
cfm.n ipn,m 



where and gm,n are some numbers, then 



Sm{Ai...A„) = 

.(1) 



2'^a!(-l)^ 



E 



(2-)! .....,.2„ 

■ 9m2a~l,m2a^rai,vi'^m2,V2 ■ ■ ■^m2a,i^2a (2''') 

This identity can be proven using Eqs. (|24II - H27|) and 

= 2'^(7!(-l)^^. (28) 

If conditions hold then again from Eqs. ll^- |[?7|) 
one can derive the identity 

B, Sm{Ai . . . i,)] = Sm{[B, ii . . . A,]), (29) 

where [. . .] denote a commutator. 

V. CLUSTER STATES AND GROUND STATE 
ENERGIES OF 

Let us consider the states 

|ncOo = Q+"|0), n = 0,...,(2? + l), (30) 

where the operator Q is given by 

Q+ = Sm{Qt^^)- (31) 

creates a (2s -I- 1) particle state. In the following 
we shall prove that \nd)^ is an eigenstate of the model 
hamiltonian hm It is easy to show that the state 

|nc^)o is annihilitated by the Casimir operator Co (See 
Eq. (HHl) oiSO[2l + l) 



Co\ncl)a = 



(32) 



by moving the operators Bl^m towards to the vacuum 
using the identity 



Bl,M, SmiQofi) 



\L-M 



2o-( 1)^ _c //)+ 



(33) 
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where the operators 



Q 



+M 



\ TO Ad — m —M I 



create also pair states with S = Q, but with angular mo- 
mentum quantum numbers L and M . For L — M = Q 
this expression agrees with (jJIJ. It is easy to prove that 

Q2p+i,M = 0, p: integer. 

Co is build up from Bl,m with L odd but the right hand 
side of in that case zero, because Q\ m = for L 

odd. As a result Co can be moved towards to the vacuum, 
which is annihilated by Co- As a side result we have the 
property 



= 0, 



(35) 



i.e., the operator Q"*" commutes with Co- 
in order to show that |ncZ)o is also an eigenvector of 
Cu it requires more elaborate calculations. In moving 
the operators B]^ m towards the vacuum one encounters 
several new objects from the commutators. For them one 
can use the identities (See Appendix 



(36) 



L=0 M=-L 



and 

21 



= -fI^™(g+,S)'- (37) 



L=0 M= 



As a result one has 



Cu\nd)o 



[I] 



m] + mnd)o, (38) 



i.e., \nd)o as given by (|Sn|l is really an eigenvector of 
Correspondingly for hm using (|32|l one obtains: 



hm\nd)o = e„|nd)o. 



(39) 



with 



n(2Z + l)(2s + 1) 



[n{2s - 1) + 2/ + 1 + 4s] . 



8(2/ + 3) 

This expression can be cast into the form 

n(n — 1) 



nei 



-1, 



(40) 
(41) 



where 



and 



(2/ + l)(2s + l) 
4(2/ + 3) 



(3s + /) 



^_(2[+l)(2s + l) 
' 4(2/ + 3) ^ ' 



(42) 



(43) 



The first term on the right hand side of Eq. H41I) can be 
interpreted as the energy of n independent one-clusters 
and the second term as a kind of cluster-cluster interac- 
tion energy. Furthermore the energy ei can be rewritten 
as 



(2^ cr((T - 1) 



where 



.(2) 



2/ -1-1 



5 = 



12 



.(2) 



2/4-3 



(44) 



(45) 



are independent of the spin. For spin one-half particles 
7 = and the prefactor of 5 in H44I) is zero, which means 
that the clusters consist of pairs which are independent. 

One can raise the question what is the ratio of the 
interaction energies of two pairs within the same cluster 
and when they belong to two different clusters. This ratio 
is equal to 



7/0-2 



3(2s + l) 
2s- 1 ' 



s > 1. 



(46) 



It is remarkable that this expression is /-independent. 
This ratio is always bigger than one, monotonically de- 
creasing with increasing s. 



VI. AVERAGE ENERGIES IN CLUSTER 
STATES 

Acting with the operators B\ (See Eqs. 0,10) on 
the cluster states of using the rules (|B1|) . H29|I one 
gets three terms: 

1 ,x 4:na\L] , inaia ~ 1) , . 
Bl\nd)o = ' \nd)o + -\aL) 



[1] 

4n(n - 1)0-2 



[I] 



(47) 



where the vectors |q!l) and \(3l) are defined as follows: 



(48) 



M=-L 



and 



\h) = o+("-^) E (-1)' 



^s^iQiMr )s^iQt-MQor M- (49) 
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In the special case L = the two vectors agree with 
\ncl)o, and for i = 2, 4, . . . , 2/ it can be shown that ja^) 
and |/3i) can be decomposed as 



4[L] 
2/ + 3 



\ncl)o 



l«L>, 



L = 2,4,. 



,21 (50) 



and 



l[s]{2l + 3) 



|nd)o + |/3f), L = 2,4,...,2Z 



(51) 

(See Appendix El, where \aj^) and |/?^) are possibly 
zero vectors, but if not, they are eigenvectors of Co with 
positive eigenvalues, and are automatically orthogonal to 
\ncl)o (which is also an eigenvector of C'o but with zero 
eigenvalue) . 

Using the above results the energy of h averaged with 
\ncl)o is 



{ncl\h\ncl)o {ncl\hm\ncl) q 



{ncl\ncl)o 



{ncl\ncl)o 



(52) 



where e„ is given by according to . 

In Appendixiniin the special case s = 1/2 we prove the 
stronger statement HC7|I , namely that the ground state of 
h is that of hm for even particle numbers. This provides 
an alternative proof within our framework of Racah's re- 
sult namely, that the ground state of spin 1/2 parti- 
cles interacting by an attractive (5-function potential con- 
sists of independent pairs. 

In order to present some numerical spectra in the most 
suitable form let us introduce the energies E' defined as 



E' 



Ns{2l + 1) 



(53) 



(|53|l makes such a shift, proportional to the particle num- 
ber N, that E'{N = 0) = E'iN = {21 + l)(2s + 1)) = 0, 
i.e., the energies E' are zero at zero and complete fill- 
ing. At the cluster numbers iV„ the particle number is 
N — n{2s + 1). Making the shift as in (|53|) the ground 
state energies E' of the approximate Hamiltonian hm at 
the cluster number 7V„ are on the curve 



E'{N) = 



N{2l + l){2s- 1) 
8(2? + 3)(2s + l) 



[(2s + l)(2/ + l)-iV] 



(54) 

We show the numerically calculated ground and all ex- 
cited state energies E' as a function of the particle num- 
ber TV, for I = 3 and s = 1/2 in Fig. ^ and for / = 2 and 
s = 3/2 in Fig. 121 

The spectra are symmetric with respect to the point 
N = {2l+l){2s+l) /2. This symmetry is a consequence of 
a particle- hole transformation • Our data show that 
the lowest energy states at particle numbers -/V„ are non- 
degenerate and there is a gap to the first excited states 
at the same particle numbers. 

The regularity in the spectrum in Fig. |21 as compared 
to the spectrum of Fig. ^ presumably arises from the 




FIG. 1: The energy levels E' of the dimensionless hamilton 
operator as a function of the particle number A'' for I — 3, 
s = 1/2. Individual energy levels are denoted by crosses (+). 
The dashed line is the function I54II . 



+ * ;; 



* + + + + + + + + 



FIG. 2: The energy levels E' of the dimensionless hamilton 
operator as a function of the particle number A'' for I — 2, 
s = 3/2. Individual energy levels are denoted by crosses (-|-). 
The dashed line is the function 15411 . 



fact that Fig. 121 refers to a situation when the hamilton 
operator coincides to the model one ()15fl while this is not 
the case in Fig. 



VII. SUMMARY AND CONCLUSIONS 

In preceding sections we have analyzed in great de- 
tails the clusterization phenomenon in the one open shell 
model with attractive 6 interactions. We have con- 
structed the states \ncl)o and showed that in a certain 
subset of shell parameters I, s these states are exact, and 
in other cases they are quite good approximate ground 
states. We also showed that there exist a cluster creation 
operator by which with repeated applications to the 
vacuum we can create the (approximate or exact) ground 
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states. We expressed the cluster creation operator in 
terms of the symmetrized product of the pair creation 
operator Qq g by introducing the symbol Sm- 

It is interesting to note that the symbol Sm itself, 
which by definition symmetrizes with respect to the an- 
gular momentum indices, can be interpreted as an an- 
tisymmctrizer which act on spin indices. This follows 
from the form l|A2|l of the cluster creation operator 
Q+ = SrniQlZ) given by Eq. 

Next, let us discuss the main analytic results (|52|l . This 
equation shows that in the cases h ^ hm, i-S., when the 
model and the one open shell Hamiltonian differs, \ncl)o 
can be regarded as a variational ansatz, i.e., the ground 
state energy of the n-cluster state lies below e„ given by 
(|40|l . In practice |nc/)o is expected to provide quite a 
good approximate ground state of h as shown in Section 
mil for a number of examples. Furthermore, if we make 
the decomposition h = hm + R and if we consider hm as 
an unperturbed Hamiltonian with unperturbed ground 
state \nd)Q then Eq. H52I) shows that we do not obtain 
correction to the unperturbed energy e„ in first order, 
corrections are at least of second order, which can be 
small because the off-diagonal matrix elements can be 
small. This explains the numerical finding for / = 3, 
s = 3/2 and = 4 that the deviation of e„ from the true 
ground state energy is of the order 0.3% (see Sec. IIIIIl . 
The true ground state, in general, is not the state which 
is annihilated by the Casimir operator of SO{2l + 1), 
but there is a small admixture which mixes some other 
eigenstates of Co with positive eigenvalues. Note that it 
remains antisymmetric in the spin indices representing a 
cluster state. 

The main technical problem in deriving values for the 
approximate (or exact) cluster energies e„ was relegated 
to Appendix ^ Here we applied a powerful projec- 
tion technics in a non-orthonormal and linearly non- 
independent basis. This procedure showed that the un- 
perturbed problem for hm is non-degenerate, if we follow 
how the Casimir of SO{2l + 1) acts in the relevant vector 
space. 



APPENDIX A: PROOF OF THE IDENTITIES 
(|36|l AND (|37|l 

Let us introduce the spin dependent antisymmetric 
matrix G^^z by 



(Al) 



If, on the left hand side of Eq. (|36|l we use once again 
(EZj) with Ai = Qt,M' ^2 = Qt.-M and with A3 = . . . = 
Act = Qto and perform the standard sum over L and M 
we arrive to 



21 L 



5Z $Z ( l)^^'5'm(QL,M*9L,-Af*9o^O ^) 



L=0 M=-L 



2^a!(-l)- 



(2a)! {[l][s]y- 



E 



£1/1,..., 1/2^ '^mi.ui '^ni, 1/2 "^-mi, 1/3 ^—711,^4 



(A3) 



If now we change the order of o,^mi using anti- 

commutativity we can get two more G factors. In the 
next step if we change the second and third indices of 
the e tensor after comparison with (|A2|) we arrive to the 
operator identity Eq. 

In order to proof Eq. (|37|l we proceed as above. Chang- 
ing the order of two creation operator (and dividing both 
sides with a common factor) it is left to proof that 



- y 



^Vi,...,V2cr Gl/i .112 C^V-i.Vi ■ ■ ■ Gjyjo- _ i , 1/2<t 



<^1,--'^2ct 
«l,---.-2<T 



XCzi,.. .,22^^2^.22 62:3,2:4 . . . >-'22ct-1,22ct 



.G2 



E 



Cfl,... ,1^2^^1/1,21^1/3,1/^ . . . Gi/2^_i,i/2„ 



'^1 ^2a 

= 1' ■■r^2(7 



XC2i,...,22ct G,/2, 22^23,24 ■ ■ • G z2^-i,Z2a ■ (-^4) 

By introducing the spin dependent matrix operator 

Fu.z ^ ^ ^t'.I'2---,l'2cr G2, 1/2 G1/3 ,1^4 ... (5*1^2^-1 ,!^2<T (-^^) 

f2,---,l'2<T 

and the spin-scalar operator 



^= E 

Vl,...,V2a 



^Ui,...,li2,,Gui.li2Gu^.lii ■ • • Gi/2„„i,l/2„ {^^) 



from (|A4|) it is left to be proven that 

1 . . . 

27Tl^'" E 

h',Z — — S 

Next, we show that 



(A7) 



If we use H27(l with Ai = . . . = = QJq have 

, tT(cr-l) 

X 01^1,1^2^1^3,1/4 . . . Gi/2„_i,l/2„ 



(A2) 



Fi/,2 = 



2.S + 1 



-F. 



(A8) 



If this is true then Eqs. I|A7|1 . IIA4p and correspondingly 
are also true. 

Let us study first the v ^ z case in Eq. (|A5(I . In 
that case among the summing indices {y^, ■ ■ ■ , i^2a) of the 
antisymmetric tensor e the index z should occur. Let it 
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be the index Vi = z. This cannot be because Gz,z = 
0. It means that two different G has the same index 
z. However, in the indices V2 and J^j+i (if the index z 
occur in the first index of the second Gz,vi^i) or z/i_i 
(if the index z occur in the second index of the second 
Gui_i,z factor) the product of two G factor is symmetric, 
while the antisymmetric tensor e in the same indices is 
antisymmetric. Summing over V2 and ViJ^i or Vi-i we get 
zero. It means, that 



F^^z =0, if 7^ z. 



(A9) 



Next, let us study J^^^ with v fixed. This particular 
index v should occur among the summing indices of F . 
It can be: vi = v or V2 = v , . . . , or V2iy — v. Changing 
the order of indices in all cases such that v be the first 
in the e tensor and at the same time in the first G factor 
(Using antisymmetry of the G-s and e) all the terms give 



F=(2s + 1) ^--^ 



J Gi,_u^ G 



Gu 



= (2s + l)F,, 



(AlO) 



Eqs. ljA9(l and (|A10|I together gives the operator identity 
which completes the proof of Eq. H37|l . 



APPENDIX B: DECOMPOSITION OF |qi,) AND 

Let us consider how la^) or behaves on applying 
Co to these vectors. Here we treat the calculation for 
\I3l)- Exactly the same method can be apphed for |q;l) 
with one minor difference, which will be shown below. 

It is easy to show that 



21 M 



B 



Li,Afi, Af2 



L=o M=~L 



[L]{-1 



M I Li L2 L 
\ Ml M2 -M 



Li L2 L 
I I I 



Q 



(Bl) 



with L2 even. In using Eqs. and HBljl one obtains 
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Co\f3L)^2[l]{l~SL,o)\l3L)-m E 



ti=o 




(B2) 



with L even. In other words, the vector space V spanned 
by the vectors \L = 0,2, . . . , 21} is invariant under 
Go- If we define matrix elements in V for an operator O, 
for which V is an invariant vector space by 



(B3) 



where X^l' stands for X]l'=o L'-even^ matrix elements 
of Co can be read off (IB2II : 



(B4) 



Co\I3l)=Y!^l',l\Pl') 



GL',L = 2m^6Lfi)SL',L~m 




L I I 
L' I I 



X 2[1]5l',l - 4[L] 



= (1-<5l,o)(1-<5l',o) 

^ 1_ / L II 
l] \ L' I I 



for L,L' even. The second equality follows from 
I M I I L M I W 



(B5) 



(B6) 



for L even. Actually the same matrix elements occur in 
vector space V = {\aL)\L — 0,2, . . . , 21}: 



(B7) 



Unfortunately, neither the vectors la^), nor the vectors 
\(3l) for L — 0,2, ... ,21 are linearly independent. This 
clearly follows from Eqs. H36() and H37I) if we apply both 
sides to the vacuum and multiply from the left by an 
appropriate power of . However, even in the linearly 
not independent (and correspondingly not orthonormal) 
basis one can still calculate matrix elements of operator 
products such as 



(B8) 



provided V is an invariant vector space of O^^-* and &'^\ 
and furthermore, the matrix elements for O^^^ and O^^^ 
are fixed (This statement follows from (|B3l) if O = 0(2) 
and we act on both sides with O^^^ from the left). If one 
calculates matrix elements of the operator R 

R = C^~ {81 + 6)C^ + 8l{2l + 3)Co (B9) 

by the well-known properties 0| of the Wigner 6j- 
symbols it turns out that 

Rl',l = 0. (BIO) 

It also means that the operator Co on V or V' fulfills 

Co(Co-Ali) (Co-{4l + 6)i) =0, (Bll) 



where / is the identity operator. 

If an operator O has the the minimal polynomial 



= (6 - Ai/) . . . (o - xj) 



(B12) 
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with finite d then O admits the decomposition 

d 

= ^A,A, (B13) 



where P, is a projector, i.e., PiPj — Pi5i,j and 



n 



3Vi 



(B14) 



From Eq. HB11|I is clear that on V at most we have 
d = 3, and the three eigenvalues of Co are Ai = 0, A2 = Al 
and A3 = (4/+6) respectively. Our main purpose is to cal- 
culate the orthogonal projection of toQ+"|0) = |/3o)- 
This latter vector is an eigenvector of Co with eigenvalue 
0, thus let us consider the projector of Pi. 



Pi = 



Co - 'ill) [Co - (4Z + 6)/ 
4/(4/ + 6) 



(B15) 



Taking matrix elements on both sides is easy. Straight- 
forward calculation leads to 



(1-<5l,o)(1-^l',o) 

(B16) 



If we consider the decomposition 

= A|/5l> + (/- A)|/3l) 



(B17) 



the first term Pi\(3l) is vector which is an eigenvector of 
Co with eigenvalue zero, the second term, if it is nonzero, 
however belong to the subspace in which Co has positive 
eigenvalues. Thus, the to vectors on the right hand side 
of l|B17p arc orthogonal to each other. Most important 
is the first term. Knowing the matrix elements of Pi it 
reads as 



(BI8) 



Equation H37|l implies (by applying both side to the vac- 
uum and multiplying by Q^'-""'^-' from the left) 



21 

E 

L' = 2 



+ M 



\Po) 



(B19) 



Putting Eqs. (|B17I) - (|B19|) together we obtain the antici- 
pated result of Eq. 

To prove the corresponding result (|50|l for the decom- 
position of Iol) we can proceed as above, but instead of 
(IB 191) we should use 



21 

E 



\aL' 



2l\ao) 



(B20) 



which follows from the operator identity (|36|l if we apply 
both sides to the vacuum and multiply by Q^'"^^^ from 
the left. 



APPENDIX C: THE SPECIAL CASE s = 1/2 

In the special case s — 1/2 the parameter a is equal to 
1 and because of the only two possibilities for = 1/2 =t 
or = — 1/2 =1 there is no need to symmetrize as 
defined in 



Consequently, according to 

By definition Eq. (jSU gives for s = 1/2 



(CI) 
(C2) 



I I L 



1L,M 
I 

E 

rn— — l 



m M — m ~M 



(C3) 



Acting with h to \nd)o Eqs. O, ^ and ^ for s = 1/2 
leads to 



h\nd)o = e„|nc/)o - —n 



21 / 

(--i)E 

L=0 \ 



/ / L 




(C4) 

Next we show that the second term is identically zero for 
s = 1/2. Using ljC2|l ljC3|l and Eq. lO it is easy to check 
that 



L=0 



+m2 



X /mi 



m,2 ;m3 ,7714 



a+ ,a+ ,a+ ,a+ ,10) (C5) 

mi,T m2,l — "13,1 -"i4,il ' \^^J 



If we insert into this equation the second form © for / 
we arrive at the rather long form 



Efo J o^\^^) = ^Qot'' E (-in 



+m3 



Z,:cvcn 

The operator t^ma i^rn^ i antisymmetric in 
TOi,m3, while its coefficient is symmetric in the same 
indices, thus, the right hand side is zero, this completes 
the proof that the state |nc/)o is an eigenstate of h: 



with 



h\ncl)Q = e„|7ic/)o, s 1/2, (C7) 
e„ = -^^, . = 1/2. (C8) 
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